Expansion of a Bose-Einstein Condensate in an atomic waveguide 
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The expansion of a Bose-Einstein condensate in an atomic waveguide is analyzed. We study 
different regimes of expansion, and identify a transient regime between one-dimensional and three- 
dimensional dynamics, in which the properties of the condensate and its further expansion can be 
well explained by reducing the transversal dynamics to a two-level system. The relevance of this 
regime in current experiments is discussed. 



During the last years, the Bose-Einstein Condensa- 
tion (BEC) H has constituted one of the most active 
research fields in modern Atomic Physics. Several exper- 
iments have shown the coherent character of the outcou- 
pled bosons from a condensate Q], i.e. a BEC constitutes 
a bright coherent matter-wave source, offering promising 
possibilities in the context of matter-wave optics. How- 
ever, contrary to the extremely dilute non-condensed 
cold atomic clouds traditionally used in Atom Optics || , 
the atom-atom interaction induces the BEC dynamics to 
be intrinsically nonlinear. In this sense, the BEC optics 
has also been called Nonlinear Atom Optics (NLAO). Re- 
markable experiments have investigated the resemblance 
between nonlinear optics phenomena and those in NLAO, 
as four-wave mixing Q and dark-solitons The non- 
linearity may deeply alter the dynamics of the wavefunc- 
tion, and the straightforward validity of atom optics el- 
ements should be questioned. As a consequence, BEC 
optics should not be considered only as a trivial exten- 
sion of atom optics, but an emerging discipline of its own. 

The analysis of low-dimensional BEC has been re- 
cently a subject of active investigation In these 
systems one or two spatial directions do not contribute to 
the dynamics, since the mean-field energy is smaller than 
the typical trapping energy in these directions. New ef- 
fects have been predicted in low-dimensional BEC, as 
quasicondensates [|l2|,|l3| (for a recent experiment, see 
Ref. ||), and Tonks gas p3| , pT[ . In addition, the physics 
of low-dimensional Bose gases is closely related with one 
of the most actual challenges of the modern Atom optics, 
namely the developing of atom optics micro-structures 
Jl5],[l6]]. In these devices, an integrated element should be 
designed to solve the fundamental problem of wiring, i.e. 
to allow for signal propagation between different parts 
of the structure, with negligible leaks. In this direction, 
current-carrying structures |l6| or dipole-force confine- 
ment in light fields ||[l5 17 constitute the most promis- 
ing possibilities. Up to now, the majority of experiments 
with atomic waveguides have used cold non-condensed 
atomic clouds. However, recent experiments in which a 
BEC have been confined into an atomic waveguide JTt| ] 
open the interesting and timely question of how a BEC 



expands inside of an atomic waveguide. It is the aim 
of this paper, to characterize in detail the loading of a 
BEC into a waveguide by the expansion of an initially 
3D trapped BEC @. The BEC properties (both the 
initial chemical potential and the subsequent expansion 
dynamics) strongly depend on the dimensionality of the 
guiding, ranging from a purely 3D character to a fully ID 
one. We shall show that due to the cylindrical symmetry 
of the guiding potential, there is an intermediate regime, 
significant for a large number of experimental situations, 
in which the BEC is well described by considering the 
radial degree of freedom as a two-level system, signifi- 
cantly departing from the ID predictions, even when the 
excited radial states are scarcely populated. As a con- 
sequence, very large radial frequencies are necessary to 
neglect observable 3D effects on the dynamics, even for 
a small number of condensed particles. 

In the following, we consider the guiding geometry cur- 
rently used in experiments performed at the University 
of Hannover |17|, although our study may be equally 
extended to similar geometries, as quantum wires [ |l6| . 
Typically, a BEC is initially created in a hybrid trap, pro- 
duced by the overlapping of a cylindrical magnetic trap, 
and a blue-detuned TEM i (doughnut-mode) laser beam 
propagating along the axis of the magnetic trap. The cor- 
responding potential exerted on the atoms can be written 
as V(y) — Vi(r) — mw 2 r 2 /2 + muj 2 z 2 /2, where uj r and 
uj z are the radial and axial trapping frequency respec- 
tively. After the creation (at t = 0), the magnetic trap 
is switched off. The subsequent evolution corresponds, 
therefore, to a BEC confined radially by the laser field 
but free to move in the axial dimension, as in a waveg- 
uide. We shall assume that the radial frequency is kept 
unchanged after the releasing of the magnetic trap, and 
therefore for t > 0, V(r) = V f (r) = muj 2 r 2 /2. To facili- 
tate the discussion of our results, we shall take as a ref- 
erence case a 87 Rb BEC of 10 4 atoms initially located in 
an hybrid trap of ui r — 2ir x 450Hz and ui z — 2ir x 10Hz. 
This is a typical situation in the ongoing experiments 
performed at the University of Hannover. 

For sufficiently low temperature, the BEC dynamics 
can be accurately described by the corresponding Gross- 
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Pitaevskii equation (GPE) m 



^ 2 + V + g\^v,t)\ 2 \r. (1) 
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where the potential term V(r, t) is defined above, and g = 
AirTi 2 aN/m is the coupling constant, with a the s-wave 
scattering length (5.8nm for 87 Kb), and m the atomic 
mass, while iV is the number of condensed atoms. With 
this choice, ip(r,t) is normalized to 1. 

The BEC expansion in the waveguide after releasing 
the axial trap is a particular case of an anisotropic har- 
monic trap with time-dependent frequencies. For suffi- 
ciently shallow traps, the BEC can be considered at any 
time to possess a self-similar 3D Thomas-Fermi (TF) 
wavefunction p9|,pfj|: 



iP(r,t) = b^ 1 b- 1 / 2 X (r/b r ,z/b z )e- lf " TW/h e-^ r ' z ' t \ (2) 

where 4>(r,z,t) = (m/2h)(r 2 b r /b r + z 2 b z /b z ), x(r, z) 2 — 
(jii — Vi(r))/g, b r (t) and b z {t) are scaling factors which 
evolve in time according to b r (t)+L0 2 b r (t) — lu 2 /bf(t)b z (t) 
and b z (t) = co 2 /b 2 r (t)b 2 z (t), and r(t) = /' dt' /b 2 r (t')b z (t'). 
We shall define the /oose-guiding (LG) regime the situa- 
tions in which the 3D TF approximation is valid. In this 
regime, the initial BEC in the hybrid trap is character- 
ized by the corresponding chemical potential 
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The TF approximation, which neglects the kinetic en- 
ergy term, is roughly correct for trapping energies smaller 
than the corresponding self-interaction. Since the trap- 
ping geometry used in the guiding experiments is strongly 
cigar-shaped, the TF approximation may become invalid 
in the transversal directions, while still being accurate 
axially. We shall refer to this situations as tight-guiding 
regime. An extreme case of this regime occurs when the 
radial binding is strong enough to dominate completely 
the transversal dynamics, which can therefore be consid- 
ered as frozen. In this case, the dynamics becomes ID 
in the axial direction, and the wavefunction can be fac- 
torized as ip(r,t) — ipiD(z,t)&o(r), where 3>o( r ) is the 
ground-state wavefunction of the transversal harmonic 
oscillator. Provided the TF approximation is valid in 
the axial coordinate, the dynamics of the BEC in the ID 
guiding limit can be accurately described by a self-similar 
TF solution in the axial direction, 



(z,t) = 



where Xid(z) 2 = {hid - V(z))/g 1D , 
mb z z 2 /2hb z , b z {t) = u) 2 z /b 2 z {t), r{t) = f dt' /b z (t'), and 
V(z) = muj 2 z 2 /2. After integrating over the transversal 
direction, one can find the effective nonlinear parameter 
9iD = g J rdrd^oir)^ = (muj r /2Trh)g. By properly 
normalizing, the ID chemical potential becomes 
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However, in a general tight-guiding situation, the 
transversal dynamics cannot be considered as frozen, 
even though the TF condition is violated in the radial 
direction. Since tvjj r > g\ip\ 2 , only a few number of ra- 
dial excited states is occupied. In particular, for a cylin- 
drically symmetric trap only excited states of zero po- 
lar angular momentum are possible, i.e. only one ex- 
cited state $„ is possible within the 2n-th radial energy 
shell. This produces that for a large number of cases of 
tight-guiding waveguides (even for \i ~ 2%uj r ), only the 
first excited state $1 is significantly populated before the 
opening of the hybrid trap. In that case the transversal 
degree of freedom can be well described by a two-level 
system formed by the first two radial eigenstates: 



^(r,t) =CoCM)$o + CiCM)*i 



(0) 



Due to the tight-guiding, |Ci|/|Cb| -C 1. Expanding up 
to first order in this small parameter, and assuming TF 
axial profile, the GPE at t = can be cast to the form 



fltC'o = V(z)Co + gmCl - ^C 2 C U 



(7) 



JhCi = {2huj r + V(z))Ci - 9 -^Cl + ^C 2 ^. (8) 

Cq and Ci are taken as real functions, and fl t = — huj r , 
with /i t is the chemical potential in the tight-guiding 
regime. Keeping only first order terms in |Ci|/|Co|, we 
find that the axial profile n{z) — J dxdy\ip(x, y, z)\ 2 can 
be written in the form 



n(z) ~ n 1D (z) 



4 \2huj r /g 1D - n 1D (z) 



(9) 



where nm {%) = |xir>(^)| 2 is the the axial density in the 
ID limit, while the second term within the brackets in- 
troduces the correction coming from the transversal dy- 
namics. The value of the chemical potential can be found 
by imposing normalization to the wavefunction: 



1 = / n(z)dz 
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where [3 — 2%u) r / (aid — 1< The root fx of this equation can 
be found easily by graphical methods. For fL t /2huj r -C 1 
one can find the analytic correction to the ID solution: 



f- — fJ-iD (1 - HlD/^TwJr) . 



(11) 



In order to test the above discussion, we have numeri- 
cally solved the 3D GPE by means of a Crank-Nicholson 
method. After obtaining the ground-state of the hybrid 
trap using imaginary time evolution, we have calculated 



2 



the corresponding chemical potential (i and compared it 
with the analytical predictions for the different guiding 
regimes (Eqs. ([§, (|) and (|l0|)), for different trap ping pa- 
rameters and number of condensed atoms. Fig. |lj shows 
the dependence of /i with respect to the radial trap fre- 
quencies, for the reference experiment discussed above. 
As expected, (||) becomes inaccurate in the tight -guiding 
regime, and vice versa (|^) provides wrong results for the 
LG region (shaded region). It becomes clear from the 
figures that the corrected chemical potential, which ac- 
counts for the two-level transversal dynamics, describes 
very well the results obtained by direct numerical integra- 
tion of the 3D GPE out of the LG region. In particular, 
it accounts very well for the transition region between 
the 3D and the ID regime, and it will merge with the 
ID result for low chemical potentials. It is particularly 
interesting that a radial compression does not easily lead 
to a fully ID behavior (even for 10 4 atoms), unless very 
large compressions (oj r > 2tt x 10 4 Hz) are applied. On 
the contrary, there is a large regime of trap aspect ra- 
tios (u> r = 2tt x 500 Hz^ 2tt x 10 4 Hz) for which Eq. 
( |l0| ) describes the actual chemical potential of the BEC. 
We should note that the ID prediction can significantly 
overestimate the chemical potential. As an example, for 
u) r = 2ir x 1.35kHz, fim is 10% larger than the actual 
one, even taking into account that more than 95% of the 
atoms are in the radial ground state. This difference will 
have significant consequences in the BEC expansion after 
the releasing of the axial trap, as discussed below. 
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FIG. 1. Chemical potential as a function of the radial fre- 
quency, u) r /2n. The LG region is shaded. 

When the magnetic trap is released, the BEC ex- 
pands axially within the waveguide produced by the 
doughnut laser, until reaching a stationary mean ve- 
locity. We have numerically simulated this process by 
solving equation (Q) in real time. We shall concentrate 
in the dynamics of the BEC variance, defined as a = 
\J (z 2 ). For the LG regime: oi (t)=ai(0)b z (t)/V7 7 where 
ct;(0) = a z (l5N(a/a z )(uj r /uj z ) 2 ) 1/5 , with a z = ^Ti/muj Zl 
whereas for the ID regime: a\D{t) = <7iD(0)b z (t)/\/E, 
with (Jid(0) = a z (3N(a/a z )(uj r /uj z )) 1 ^ 3 . Therefore, at a 
given time, v = da z /dt has different laws for the differ- 
ent regimes: vi ~ N 1 / 5 ^^ and v\u ~ N 1 ^^ 3 ■ Fig. ^ 
shows the variation of v at 50 ms, for the cases considered 
in Fig. [I]. At this t, v has reached an almost stationary 
value. From the figure, it becomes clear the transition be- 



tween the 3D behavior and the ID one. However, as in 
Fig. [I], the regime for which the 3D effects on the dynam- 
ics become unobservable is not easily achievable by in- 
creasing uj r . This point becomes clear in the inset of Fig. 
||, where the numerical results are shown to be parallel to 
the ID ones for a large range of radial trap frequencies. 
In particular, v can be significantly lower than v\d even 
for large compressions. From the results of the previous 
section, this effect can be explained by (i) the reduction 
of the initial chemical potential due to the non-negligible 
influence of the population of the first transversal state, 
and (ii) the departure from the TF profile due to the 
same reason. As an example, for oj r = 2tt x 1.35kHz, at 
50 ms vid — v ~ 0.15 mm/s, which determines that the 
condensate variance becomes approximately 7/im smaller 
than that expected from a ID solution. Such difference is 
experimentally observable with the current imaging tech- 
niques. In the two-level regime, Eq. (|l|) is equivalent to 
the system of time-dependent differential equations: 

- gMC^Cx - \gi D ClCt, (12) 



■^9id\Co\ 2 C + -gi D ClC{. 



(13) 



The numerical simulation of the previous equations, also 
depicted in Fig. ^, shows an excellent agreement with the 
direct integration of Eq. (|l|). Therefore, the numerically 
hard task of calculating the 3D corrections by solving the 
3D GPE, can be performed in a much simpler way, by 
simply solving a set of 4 ID coupled differential equations 
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FIG. 2. Velocity of expansion after 50ms as a function of 
the radial frequency, u) r /2n. The LG region is shaded. Inset: 
detail of the tight -guiding region. Hollow squares indicate 
the solution of (^) and (jl|) in the the two-level regime. 

Finally we shall analyze the acceleration of the expan- 
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sion. It is easy to find that for the ID regime 
3 hu r aN 1 
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On the other hand, since u r > lu z in a general guiding 
situation, we may consider b r to follow adiabatically b z , 
and in the LG regime b r (t) — 6 z (i) _1 ^ 4 . Then, 
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(15) 



Therefore, in general, the acceleration decreases as Hit) ~ 
l/<j n (t), with n ranging from n = 3/2 in the LG condition 
and n — 2 in the ID regime. Fig. ^ shows n obtained from 
the 3D GPE, for the cases of Fig. 
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FIG. 3. Coefficient n (see text) as a function of the radial 
frequency, uj r /2n. The LG region is shaded. 

In the present paper, we have analyzed the chemical 
potential, and characterized the expansion dynamics of 
a BEC within a waveguide. We have numerically in- 
vestigated the transition between a 3D situation (loose 
guiding) and a ID one. Although the extreme 3D and 
ID situations are well described by the corresponding 
TF self-similar solutions, there is a broad region of in- 
termediate parameters which is not well accounted by 
these solutions. Due to the cylindrical symmetry of the 
guiding potential, we have shown that in this region the 
transversal dynamics of the BEC can be modeled as a 
two-level system, formed by the ground and first excited 
state of the radial trapping potential. This model is in 
good agreement with the direct numerical simulation of 
the corresponding 3D GPE. As we demonstrate, even a 
small population in the first radial excited state can in- 
troduce observable differences between the actual situa- 
tion and the ID one, in the resulting chemical potential, 
and in the subsequent expansion of the BEC after releas- 
ing of the axial trap. From our results, we can predict 
that a very tight transverse confinement is needed to ne- 
glect observable 3D effects determined by the transversal 
two-level system. In particular, this two-level structure 
could have important effects in the reflection and beam- 
splitting of BEC in waveguides, due to the enhancement 
of the nonlinearity at the edge of the optical elements, 
which could re-couple the transversal levels. The prop- 
erties of a BEC under such conditions will be the subject 
of further investigations. 
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